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Abstract 

In this paper, we investigate some of the classical restricted three body problems at 
the nanoscale, such as the circular planar restricted problem for three Cgo fullerenes, and 
a carbon atom and two Cgo fullerenes. We model the van der Waals forces between the 
fullerenes by the Lennard-Jones potential. In particular, the pairwise potential energies 
between the carbon atoms on the fullerenes are approximated by the continuous approach, 
so that the total molecular energy between two fullerenes can be determined analytically. 
Since we assume that such interactions between the molecules occur at sufficiently large 
distance, the classical three body problems analysis is legitimate to determine the collective 
angular velocity of the two and three Ceo fullerenes at the nanoscale. We find that the 
maximum angular frequency of the two and three fullerenes systems reach the terahertz 
range and we determine the stationary points and the points which have maximum velocity 
for the carbon atom for the carbon atom and the two fullerenes system. 
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1 Introduction 



In 1966, Jones [la] suggested in the New Scientist under the name of Daedalus that the 
hollow all-carbon cage molecules were possible. However, it was not until two decades 
later that Kroto [3] experimentally discovered Ceo and C70 fullerenes by analyzing the 
resulting mass spectrometry. The dominant peak represents Ceo fullerene followed by 
a second peak representing C70 fullerene. In particular, Ceo fullerene has Ih symmetry, 
which has a spheroidal shape. 

Various nanoscale gigahertz oscillators have been proposed in the recent literature 
including the oscillation of a Ceo fullerene or a closed carbon nanotube inside an open 
ended carbon nanotube. Such oscillators have also been confirmed by experiments on 
multi- walled carbon nanotubes. Cumings and Zettl 0] remove the cap from one end 
of the outer shell and attach a moveable nanomanipulator to the core of a multi-walled 
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nanotube through a high-resolution transmission electron microscope. They observe that 
the extruded core, after release, quickly and fully retracts inside the outer shell due to 
the restoring force resulting from the van der Waals interactions acting on the extruded 
core. This experimental result lead to the molecular dynamics studies of Zheng and 
Jiang 29|] who show that the oscillating of the inner shell between the open ends of the 
outer shell of a multi-walled carbon nanotube generates a gigahertz frequency. Molecular 
dynamics simulations on gi gah ertz oscillators have also been undertaken by Legoas et al. 



19(1 and Rivera et al. 2^, |25|. In terms of mathematical modelling, Baowan and Hill 



l| investigate the force distribution for a double-walled carbon nanotube oscillator by 
adopting the continuous approach for the Lennard-Jones potential together with Newton's 
second law, assuming a frictionless environment, to investigate the associated mechanics. 
They obtain an analytical expression for the interaction force and their model also predicts 
the gigahertz oscillatory behavior for the double-walled carbon nanotube oscillators. 

Based on the molecular dynamics simulations of various nano gigahertz oscillators, Cox 
et al. 0, 0] develop an accurate mathematical model employing fundamental mechani- 
cal principles and classical applied mathematical techniques to determine an acceptance 
condition and the suction energies of a Ceo fullerene entering a nanotube. They then 
determine the minimum radius of a carbon nanotube for which the Ceo fullerene will be 
accepted from rest and the maximum total kinetic energy once the Ceo molecule is sucked 
inside the nanotube by the van der Waals forces. In addition. Cox et al. [1, 0] show 
that the gigahertz oscillatory behavior arises from the two peak-like forces operating at 
the nanotube's open ends. The analytical model of Cox et al. [2, is also extended 
to examine some more complicated structures of the gigahertz oscillators, including the 
nanotube bundle oscillators, for which a single nanotube or a Ceo fullerene oscillates inside 
the bundle 0, 0]. For other types of nanoscale oscillators, Hilder and Hill [l^, 14] find 
that the gigahertz frequencies can also be obtained from a sector of a nanotube orbiting 
inside a carbon nanotorus and an atom and a Ceo fullerene orbiting inside a nanotorus. 

Chan et al. |2|] have investigated various two body problem, namely fullerene-fullerene, 
fuller ene-carb on nanotube at the nanoscale. In this paper, we investigate some of the 
classical restricted three body problems at the nanoscale. In particular, we consider 
three Ceo fullerenes interacting with each other and initially located at the vertices of 
an equilateral triangle. The van der Waals interaction energy is modelled using the 6-12 
Lennard-Jones potential and the continuous approach for which we assume a uniform 
distribution of carbon atoms on the surfaces of the carbon nanotube and the fullerenes. 
We determine the motion of the interacting molecules as they move relatively with respect 
to each other under the influence of their mutual centripetal force. For the case under 
consideration, the circular orbiting frequencies reach the gigahertz range. We comment 
that this paper ignores any thermal fluctuations arising from the environment. 

In section 2, we develop the mathematical model for the restricted three body problems 
at the nanoscale. In section 3, we discuss the restricted problem for the three fullerenes 
system. In section 4, we investigate the capture of a carbon atom by two fullerenes and 
we present some conclusions in the flnal section. 



2 Theory 

In this section, we derive the theoretical basis for the study of the three body problems, 
namely that of three Ceo fullerenes at the nanoscale (see Fig. [T] for details). Since we only 
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Figure 1: Three body problems at the nanoscale. Fullerene centers initially lying at the vertices 
of an equilateral triange 



consider the residual molecular interactions between molecules, it is valid to examine the 
classical three body problems to investigate the same problem at the nanoscale [s, 10, 12, 



21 



271]. We model the residual molecular interactions by van der Waals forces, 



and in particular the Lennard- Jones potential [16], which is given by 



V{p) = 4e 



12 
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P 



12 • 



(1) 



where p, e and a denote the distance between two arbitrary carbon atoms, the potential 
well depth of two carbon atoms and the parameter that is determined by the equilibrium 
distance respectively. In addition, A and B denote the attractive and the repulsive con- 
stants respectively. We assume that the carbon atoms are evenly distributed over the 
surfaces of the molecules such that the pairwise molecular energy between two fullerenes, 
i.e Ylik Ylij ^(.Pkj) can be approximated by the continuous approach, which has been pre- 
viously used by Cox et al. [1, [H]. That is. 



A 



where nf,p, dSi and dS2 denote the atomic number density, the distance between the cen- 
ters of two fullerenes, the surface area element of the two fullerenes respectively. Following 
the same calculation as that derived by Cox et a/[l, [H], we obtain 



V{r) = -QQir) + Qu{r), 



(3) 
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where Qn is defined by 



Qn{r) 



r(n-2)(n-3) U2a + r)"-3 (2a - r)'^-3 r'^-S 

where a, Cg and C12 denote the radius of the Cgo fullerene, A and B respectively. Given 
that, the mutual force, F{r) between two fullerenes can be determined by 

Fir) ^ (4) 

dr 

where d/dr denotes the derivative in the r-direction. 

According to Newton's second law, we can determine the forces acting on a fullerene 
by the other two fullerenes as 

3 — 

mkfl = ^ F{rkj)— , (5) 

where m^, r^, r^j denote the total mass of the A;- fullerene, the position vector of the 
/c-fullerene and the relative displacement of any two fullerenes k and j. We note that 
Eqs. ([5]) can not be generally solved analytically. If we make f = — fl and we assume 
that R denotes the center of mass of mi and m2 allocated at the origin of the coordinate 
system. We further define p = — R = Mfi~^r^, where fi = mi + m2- Under these 
assumptions, we can write these relative displacements as 

- r"! = f , r5 - n = p + m2/U~ V, rs - = mi/x~V. (6) 
The equations of motion in terms of the Jacobi coordinates, r and p are then given by 

F{n2) ^ , 1 F(r23) ,^ 1 nnsl 

r -\ [p — mip r) [p + mip r), 



r = — 

mim2 ri2 m2 ' mi rn 



P = fJ- 

ma 



-{p + m2P r) -\ (p-mip f) 



(7) 



^'is r23 

3 Three identical Ceo fullerenes 

We now make the following three assumptions, referred to as the circular planar restricted 
problem, to determine the motion analytically: 

1. The three fullerenes are moving uniformly in circular orbits 

2. They are orbiting collectively in the same plane 

3. Their angular velocities are the same 

Upon assuming the three fullerenes are moving in plane, we can denote {xi,yi,0) by 
the coordinates of the z-fullerene. Then, the second derivative of Xi and yi can be written 
as 

Xk = mj^ l^— [Xj-Xk), yk=m^ 1^— iVj-yk)- (8) 

. ' jfc . ' jk 
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We now assume that the three fullerenes have the same (cohective) angular velocity oj, 
and they are moving uniformly in a circle, such that we can write down the rotating 
coordinate system 

Xk = £,kCOSUjt — rjksmojt. (9) 

We note that can be defined in a similar way as x^. We then substitute Eq. ([9]) into 
Eqs. (I8j)i to obtain 



ik - 2ujrjk - u?ik = "ifc ^ X] ~ 



j^k -"^ 

rjk + 2ujik - uj'^Vk = m'^ ^ —{Vj - %)• (10) 

j^k 

We can then recast Eq. (fTOl) into a single equation by using Zk = £,k + ir]k, where i denotes 
the usual imaginary unit \/— T, so that, 

z'k + i2ujZk - oJ^Zk = m^^ ^ —{zj - Zk). (11) 

. 1~jk 

Since the fullerenes appear to be at rest in the rotating frame, we can simplify Eq. (|lip 
to become 

-Zk = Xky2^^izj-Zk), (12) 

where Xk = l/{mkOj'^). If we let [5j = Xk{F{rme)/i"me), where m,£ j, then Eq. (fT2]) gives 
the following a system of three linear equations 



' {1-P3-P2)Z1+P3Z2 + P2Z3 = 0, 
< P2Zl+PlZ2 + il-f32-Pl)z3 = 0, (13) 

mizi + m2Z2 + mszs = 0, 

For the three fullerenes moving in the same plane, geometrically, two possibilities 
arise, namely; the three fullerenes are not on the same straight line or they lie on the 
same straight line. For the first scenario, the coefficients of Zj of Eqs. (llSh i . (jlSP o and 
(|13p ^ must be proportional to each other. Given that the case, we automatically have 
mi = 7712 = n^3 = nT-, which corresponds to the case of the three identical fullerenes we 
study in this paper. In addition, we require /3i = = /^s = 1/3 resulting to the three 
identical fullerenes allocated at the vertices of an equilateral triangle, i.e., r^j = r. We 
have 



mu;^ r 6 V mr 

On the other hand, we can directly obtain the same result by investigating the equation 
of motion given by Eq. (^i- Upon assuming F{ri2) = F{ri2) = F{ri3) = F{r) and 
mi = m2 = m3 = m, we deduce 
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2m F(r) ^ 1 F(r) / m m \ ^ 1 Fir) / m m \ ^ 3F(r 
_^_^f J — r-\ r — r-\ \ f = ^ 



r = 77 r H I r H I r I r H 1 I r = r, 

r m r \ 2m 2m y m r \ 2m 2m/ mr 

(15) 

which represents a simple harmonic motion with angular velocity lv, given exactly by 
Eq. (IIID2. 

Radius of Cqq a = 3.55 A 

Carbon-carbon bond length a = 1.421 A 

Mean surface density of Cgo '^z = 0.3789 

Mass of a single carbon atom rric = 1.993 x 10~^^ kg 

Mass of a single Cgo fuUerene m = 1.196 x 10"^"^ kg 

Attractive constant A = 17 A eVx A*^ 

Repulsive constant _B = 29 x 10^ eVxA^^ 

Table 1: Numerical values of constants used in the model 



Here, we carry out the theoretical results derived above. All the necessary constants 
used in this paper can be found in Table. [TJ We note that while the molecular interactions 
are attractive at long range due to the van der Waals forces, they are repulsive at short 
range as the result of overlapping electron orbitals. The total pairwise molecular inter- 
actions between two Ceo fullerenes, which is approximated by the continuous approach, 
termed the total molecular energy, is shown in Eq. ([2]) and is plotted in Fig. [2l We find 
that the most stable configuration for two Cgo fullerenes are at 9.8 A with a much higher 
energy well depth in comparison to that of the two carbon atoms that is 0.28 eV. Given the 
total molecular energy, we can determine the total molecular forces between two fullerenes 
by differentiating, and the force F(r) given in Eq. Here, we define positive forces to 
be attractive forces and vice versa. Since the centripetal forces between three fullerenes 
are provided by the attractive van der Waals forces between molecules only, upon using 
Eq. (jl4p . we can predict the collective angular velocity to, which is plotted in Fig. [H 
Since only attractive van der Waals forces contribute to the orbiting behavior, no orbit- 
ing phenomenon can be observed for r < 10 A, where the van der Waals forces become 
repulsive. Given the orbiting regimes, the maximum angular velocity is found to occur 
at 0.12 THz, which is substantially higher than all gigahertz nanooscillators proposed in 
the existing literature. As such the three Cqo fullerenes might be used as a possible ultra 
high frequency nano-device. 



4 Capture of a carbon atom by two fullerenes 

In this section, we investigate the motion of a carbon atom under the infiuence of the 
potential energy from two Geo fullerenes. We assume that the molecular interactions of 
the carbon atom is so small that it cannot infiuence the dynamics of the two Gqo fullerenes. 
From which, the Jacob coordinates, i.e. Eq. ([7]) reduce to 
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r(A) 

Figure 2: Total molecular energy V{r) of two Cgo fuUerenes 
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r(A°) 

Figure 3: Total molecular force F{r) acting on one fuUerene 
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r(A°) 

Figure 4: Collective angular velocity u versus r, where the maximum angular velocity occurs 
at 0.12 THz at approximately r = 10.5 A 



2 F(ri2) ^ 
r, 



m ri2 



rrir 



Pi 2 p2 2 



(16) 



where m, rric, pi and p2 denote the mass of the Cgo molecule, the mass of the carbon 
atom, the distance between the atom and the first fullerene and the distance between 
the atom and the second fullerene respectively. Eq. (fT6|) i constitutes the usual two body 
problem, which has the usual physical interpretation. If we assume the same scenario 
as in Section 2 but for the two fullerenes, that is the two fullerenes move collectively in 
a circular orbit of radius ii in a plane and with the same angular velocity lj, Eq. (I16p i 
becomes simple harmonic motion with angular velocity equal to ^j2F{r) jmr for a fixed r 
and this is plotted in Fig. [5l In the later part of this section, we denote R = 10.5 A to be 
the circular radii, where the collective angular velocity of the two body systems attaches 
maximum value. 

Similarly, we assume that the center of mass of the two fullerenes are located at the 
center of a cartesian coordinate system and we introduce a rotating frame (^, r/) with the 
same angular rate as the collective angular velocity w so that the two fullerenes remain 
stationary with respect to the rotating frame. The position of the two fullerenes can then 
be determined by solving + = and ^2 — ?i = simultaneously and the situation is 
shown schematically in Figl5j Since the dynamics of the two body systems is determined, 
the only problem now left is the motion of the carbon atom about the center of mass, 
which has to be solved from Eq. (jlGp ?. 

Now making the same circular planar assumptions as above, the third component, i.e 
the equations of motion of the carbon atom, i.e. Eqs. ([TO]) becomes 
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r(A°) 

Figure 5: Collective angular velocity u versus r, where the maximum angular velocity occurs 
at 0.1 THz around r = 10.5 A 



Y 



{-R/2,0) 




Figure 6: Capture of carbon atom by two identical Cgo fuUerene 
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rric [ pi P2 J 

V + = 1- |^(,, - ,3) + ^iV2 -Vs)], (IT) 

rric [ pi pi J 

where rric, Pi = y^(^ + R/2Y + and p2 = — R/'^Y + denote the mass of 
the carbon atom, ^3, r/3, ri3 and r23 respectively. For simphcity, we only work on Eq. (fT7|) i 
and Eq. (jl7p 9 can be determined in a very similar way. Starting from Eq. (jlTP i we may 
deduce 



nic \ pi 

1 fFjPi) 
rric I Pi 



R 
2" ' 



+ 



+ 



F{P2) 

P2 
F{P2) 
P2 



R 
2" 



R 
2" 



1 rav>i)9pi , g^(P2)g/02 



m, 



\ dpi dp2 dC 



1 jdVip)} 
rricX di J 



where V{p) = V{pi) + V{p 



2)- 



where <1) = $1 + $2, and <I>i : 
conclude that 



We can further simplify Eq. (llSp to become 

•■ 1 a$ 

= y(p) + mcw2^2 ^2 _ T/^.^ ^ _ ,.,2^2 



y(p) + rUcUJ r] 



1 5$ 



2mc 5^ 

and Eq. (j20p can then be recast in the form 



f/ + 2u;^ 



1 

2mc dr] ' 



dt 



1 

—rric 
2 



and on integrating both sides, we obtain the energy integral 

1 



-rric 



$ + c > 0, 



(18) 
(19) 

We therefore 
(20) 

(21) 
(22) 



where C denotes the Jacobi constant. As in the classical three body problems, we are 
interested in determining where the carbon atom is stationary with respect to the rotating 
frame under the influence of the two fullerenes. In other words, we require the velocity 
in both the ^-direction and the r/-direction to be zero, i.e. ^ = ?) = and hence from 
Eq. ([20]), we have d^/d( = d^/dr] = 0. 

The positions, for which ^ = r) = is satisfied, are numerically obtained and shown 
in Fig. [71 We note that two minima are located inside the two fullerenes, which indicate 
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Figure 7: Stationary points of carbon atom with respect to the rotating frame , 77), which are 
marked by V and the maximum velocity of the carbon atom, which are marked by A. The two 
circles represent the two Ceo fuUerenes 



the encapsulation of the carbon atom inside the fullerenes and due to the symmetry of 
the system, all minimum occur at either ,^ = or r/ = 0. Since C is an integral constant 
throughout the conservative system, we can obtain the value of C from a particular case, 
the simplest one is (^,77) = (0,0) when (■^,r/) = (0,0). C = — $(C = 0, ry = 0) is found 
to be 9xl0~^^ J. Upon substituting C back into Eq. (j22p . we can determine the velocity 
profile of the carbon atom and hence the maximum velocity of the carbon atom, which is 
also shown in Fig. [71 



5 Conclusion 

We adopt the continuous approximation together with the Lennard-Jones potential to 
model the van der Waals between fullerenes. The total pairwise molecular interactions 
of the atoms between two fullerenes, using molecular dynamics simulations usually takes 
a substantial computational effort to implement. Here, such a total molecular energy is 
computed analytically and the corresponding total molecular force can be determined by 
differentiating. For the restricted three body problems derived in section 2, we can esti- 
mate the collective angular velocity of the two and three fullerenes systems. In particular, 
the maximum angular velocity reaches the terahertz range. In addition, we determine the 
stationary points and the points with maximum velocity of the carbon atom captured by 
two oscillating fullerenes. 
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